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<D • Abstract 

The formulation of the variational problems for the solute transport 
in a fluid layer in presence of double-diffusive thermal convection is dis- 
cussed. It is shown that the variational functional obtained by Strauss 
can be generalized and the general functional leads to accurate upper 
bounds on the solute transport for the case of small and intermediate 
values of the Rayleigh number. The general functional however is a 
non-homogeneous one but for asymptotically large Rayleigh numbers 
it converges to the Strauss approximation. Thus for small and inter- 
mediate values of the Rayleigh numbers one should use the general 
functional and for vary large values of the Rayleigh numbers one can 
& ' use the functional of Strauss. 

H 

Double-diffusive comvection is of considerable interest for oceanology as it 
is connected to the convective motion of water containing some amount of 
salt. In this case of convection in addition to the temperature gradient there 
is also a gradient of the salt concentration and this complicates the model 
equations. Below we discuss the variational problem for the upper bounds 
on the solute transport in a double diffusiove convection. Our discussion will 
be based on the optimum theory of turbulence which has been developed 
in the pioneering works of Howard, Busse, Doering and Constantin 1-4 . 
We shall not discuss the Doering-Constantin approach which easily leads to 
upper bounds if appropriate background fileds are chosen 5 . Instead of this 
our attention will be concentrated on a discussion of the work of Strauss 6 
which is based on the Howard-Busse method which has many application to 
systems with thermal convection 7-11 . In order to avoid the complication of 
the analysis we shall not discuss the case of presence of rotation 12 " 16 and 
instead of this we shall follow the guidelines from 

Below we shall model the double-diffusive convection in a horizontal fluid 
layer by means of the equations of Boussinesq approximation 6 ' 19 

(1) +u-Vu = Vp - /3gSe z + jgTe z + uV 2 u 

dt po 
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(4) V • u = 

where S is the concentration field which is divides as a horizontally 
averaged part 50 El and fluctuating part S 

(5) S(x,y,z,t) = S(z) + S(x,y,z,t) 

The same division is performed for the temperature field T too 

(6) T(x,y,z,t) = T(z) + T(x,y,z,t) 

In the above equations u = (u,v,w) is the fluid velocity (which is a fluctu 
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ating quantity as we assume that there is no mean flow: U = 0[j). po is 
the average fluid density; p is the deviation of the pressure from hydrostatic 
pressure field which would exist if S = S and T = T; 7 is the coefficient 
of thermal expansion; f3 is the density change due to the unit change of 



1 In this paper we shall use two kinds of averages: 

1. Average over horizontal plane (horizontal average) 

Q(z,t) = i 1 ™ o Tj2 \ \ dxd y QixiViZ't) 

J — OC J —OO 

2. average over the volume of the layer (volume average) 

^ C co c oo pi 



{Q)(t)= lim — — / / / dxdydz Q(x,y,z,t) 
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In additition we shall assume that we investigate the process of convection under 
quasistationary condition which is defined by the requirement that the horizontal averages 
of all quantities are time independent and the horizontal averages of all fluctuation parts of 
the quantities are 0. Because of this in the text below the horizontal and volume averages 
of the studied quantities depend only on z and the horizontal averages of the fluctuating 
parts of the studied quantities are set to 0. 

3 The basis of this assumption is in the nature of the convective motion in a fluid layer. 
When the heating of the layer is not very intensive the heat is transported through the 
layer by means of thermal conduction. Thermal conduction means that there is some 
temperature distribution in the fluid layer ( i.e., T/fl and the temperature fluctuations 
are zero) but there is no motion of fluid (i.e. the mean fluid velocity U = 0) and the 
fluctuations of the fluid velocity are too. The arising of the convective motion of the fluid 
complicates the heat transfer. The heat is transferred not only by thermal conduction but 
also by thermal convection. This means that non-zero fluctuations arise and in addition 
the amplitude of these fluctuations can be large: they are nonlinear fluctuations which in 
general can not be treated by linear model equations. Because of this the obtained below 
model Euler-Lagrange integral-differential equations will be nonlinear. 
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solute concentration; g is the acceleration of the gravity; v is the kinematic 
viscosity; ks is the solute diffusivity; kt is the thermal diffusivity and e z is 
the unit vector in the direction opposite to the direction of gravity. 

Below we shall obtain several integral-differential consequences of the 
Boussinesq equations. In the proces of obtaining of these relationships we 
shall use actively the condition of incompressibility of the fluid (the conti- 
nuity equation @) and we the assumption that we study the system long 
after the last external influence of its dynamics, i.e., when the condition 
of quasistationarity is fulfilled that leads to independence of the horizontal 
averages of the flow quantitities on the times and in addition leads to the 
condition that the horizontal averages of any fluctuation quantities (and of 
combinations and derivatives of such quantities) are equal to 0. 

Taking into account the above mentioned conditions after averaging ([2]) 
with respect to horizontal plane we obtain the relationship 

,~ d _= d 2 S 

(7) Tz wS = Ks d^ 

The horizontal average of ([3]) leads to 

d^= d 2 T 

(8) d~z WT = KT ^ 

The integration of and ([8]) with respect to z leads to the relationships 

— P7 dS 

(9) wS — ks—, — = const 

dz 

(10) wT — kt~t~ = const 

dz 

The volume average of Q and (fTUj) leads to 

(11) {wS) — nsAS = const • d 

(12) (wT) - k t AT = const • d 

where AS and AT are the differences in the salt concentration and temper- 
ature difference between the top and bottom borders of the fluid layer. As 
the integration constants are the same in ([9]) and (jlip as well as in (|10p and 
(|12p we can eliminate the constants and in such a way we arrive at the two 
relationships 

-5 dS (wS) AS 

(13) W S-K S — = — K S — - 

dz d d 

(14) wT - K T — = — K T — — 

dz a d 
Below we shall perform non-dimensionalization of the quantities on the basis 
of the following units: AT will be the unit for temperature; AS will be the 
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unit for solute concentration; the layer thickness d will be the unit for length; 
d 2 1 'ks will be the unit for time and Ks/d will be the unit for velocity. 

As next step we shall ontain two relationships known also as power in- 
tegrals. Let us multiply fl2J) by 5 and average over the fluid layer. The non- 
dimensional result is 

(15) (wS) 2 - (^S 2 ) = (wS) + (| VS | 2 ) 

The second power integral is obtained after multiplication of (pQ) by u aver- 
aging the obtained result over the fluid layer and non-dimansionalization of 
what is obtained. Thus we arrive at the relationship 

(16) - R(wS) + —RaiwT) - (| Vu | 2 ) = 

where R = P g ^ Sd is the solute Rayleigh number and Ra = 7gATrf is the 
thermal Rayleigh number. 

The goal of the methodology is the obtaining of upper bounds on the 
solute transport 

(17) Nu s = 1 + (-wS) 

We shall write Nus in the form Nus = 1 + y, i.e., when Nu$ has a max- 
umum T will have a minimum. In order to come to such a form of the 
relationship for the solute transport we rescale w, S, and 6' = as fol- 
lows 

(18) w = Aw; S = BS; 6' = AO 

where A and B can be determined by substitution of (|18p in the two power 
integrals above. The result of this is 

(19) B = A Ra ^-^^=aA 
V 1 R(wS) 



(20) A 



1/2 

(wS) + a(\ VS I 2 ) 
a{(wS) 2 - (^§ 2 )). 



Now after some calculations we obtain the following relationship for T (we 
have set S = S/R and in the final result below we omit the tilde and the 
hat signs) 

(21) T = A(| V5| 2 )[(| Vn| 2 ) - Ra(w9)} + - (wS)) 2 



(wS) 2 

On the basis of (|21|) we shall formulate the variational problem below. But 
before this we shall obtain the particular case of (f2T|) discussed by Strauss 6 . 
In this particular case there is a relationship between w and 6 and because 
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of this the variational problem becomes simpler and the number of the cor- 
responding Euler - Lagrange equations decreases. In order tho obtain the 
relationship between the two above-mentioned fields we start from Eq. ([3]) 
perform a non-dimensionalization and discuss the quasi-stationary approxi- 
mation. What is obtained is 

(22) TW K + T u. VT = fl + V 2 T 

dz dz z 

where r = ks/kt- Here we shall use an approximate relationship for ^ 
which will lead to the approximation used by Strauss. We start from ([14 
perform a non-dimensionalization and obtain 

dT 1 ks 



(23) Tz = 1 - f^T) - (»r>] 

Now the assumption of Strauss is that r = ks/kt « 1 and that the 
expression in [...] of Eq. (|23p is not larger than 1 for any z. Then 

and after the approximation ru • VT m the Eq. ([22]) can be reduced to 

(25) tw = V 2 T 

and because of the fact that 9' = T/r this lead to 

(26) w = V 2 6' 

Let 9 = Ra9' and we consider the particular case of 1 — a solution of the 
variational problem 

(27) w = w(z)f(x, y); 9 = 9(z)f(x, y); V\f = -a 2 f 
Then from (|26p we obtain 

w(z)Ra 



(28) 9{z) 



a 2 



Then the relationship for T in the approximation of Strauss is 
A(| VS | 2 )[(| W | 2 ) + ^(w 2 )] + j^S) - (wS)) 2 



(29) T 



Note that the functional of Strauss (|29|) is homogeneous one which will lead 
to homogeneous Euler-Lagrange equations and will simplify their numerical 
and analytical asymptotic solutions. 
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Let us now discuss a more general approximation and its consequences 
for the variational problem, let us again assume tu ■ VT ~ and let us 
make Eq.([8]) dimensionless. The result is 

(30 » W" T = ^ 

Thus after a substitution of Eqs. (|23p and ()30p in Eq. ()22p we obtain the 
following relationship 

(31) tw{1-t[{wT) -U/T]} = t-^-'wT + \7 2 T 

Now we see that one can obtain the approximation of Strauss if one neglects 
the terms of order 0(r 2 ) in the left-hand side of Eq. (|3ip and in addition 
one has to make so-called internal layer approximation in the right-hand 
side of Eq. (|31|) . The internal layer approximation means that one has to 
assume that -^wT ~ in the fluid layer. Strictly speaking the internal layer 
approximation is valid only in the internal sub-layer of the fluid layer. U 
to some extent it is valid in the two intermediate layers of the fluid layer 
The internal layer approximation is not valid in the boundary layers of the 
fluid where all the fields (velocity, temperature and concentration) have to 
adjust their values to the boundary conditions on the borders of the fluid 
layer. When the Rayleigh numbers connected to the discussed problem have 
small and intermediate values then the thickness of the boundary layers of 
the fluid layer is large and because of this the approximation of Strauss is 
very crude. But when the values of the Rayleigh numbers become larger 
and larger then the boundary layers of the different fields become thinner 
and thinner and the approximation of Strauss can became reasonable one. 

On the basis of all this we can formulate the variational problem for the 
upper bounds on the double-diffusive convection as follows 



Case of arbitrary values of the Rayleigh numbers (Variational 

problem 1) 

Given Ra and A > 0, find the minimum M(A) of the functional: 

_ A(| VS | 2 ) [(I W | 2 ) - Ra(w9)] + (gjjg) - (wS)) 2 

(wS) 2 

Aw- -^-V 2 e - -?-w6 + ^-wUwO) - 
I /h' Ra Ra J 

4 The intermediate layers of the fluid layer are layers where the velocity field, the tem- 
perature field, and the concentration field make a transition from their almost constant 
values in the internal layer of fluid to the sharply depending on the coordinate values in 
the boundary layers of the fluid. 
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within the class of fields satisfying the boundary conditions and the conti- 
nuity equation V ■ u = 0. 

Above n is a Lagrange multiplier by means of which one takes into an ac- 
count Eq.([3TJ- 

For asymptotic large values of the Rayleigh numbers we can use the 
approximation of Strauss: 

Case of asymptotic large values of the Rayleigh numbers 
(Variational problem 2 - approximation of Strauss) 

Given Ra and A > 0, find the minimum M(A) of the functional: 

A(| VS | 2 )[(| W | 2 ) + + ((g) - (wS)f 

(wS) 2 

within the class of fields satisfying the boundary conditions and the conti- 
nuity equation V • u = 0. 

As concluding remarks we make several notes with respect to the two 
variational problems. 

• The more general variational problem 1 is a non-homogeneous one. 
Because of this the corresponding Euler-Lagrange equations will be 
non-homogeneous too which means that the task for the obtaining 
upper bounds on the solute transport will be very difficult. However 
the obtained bounds will be very accurate especially for the case of 
small values of the Rayleigh numbers where the Strauss approximation 
is relatively inaccurate because of the large thickness of the boundary 
layers of the fields of the velocity, temperature and concentration. 

• For large values of the Rayleigh numbers the upper bounds connected 
to the more general variational problem 1 will come close to the upper 
bounds of the less general variational problem 2 of Strauss and then the 
more simple and homogeneous Euler - Lagrange equation connected 
to the variational problem of Strauss can be used for calculation of the 
upper bounds. 

• The variational problem of Strauss is derived only for the simplest pos- 
sible form of the studied fields: fields with one characteristic wavenum- 
ber. The more general variational problem 1 is suitable for investiga- 
tion of bounds for optimum fields that have arbitrary number of char- 
acteristic wave-numbers (so called multi-wave-number solutions of the 
Euler - Lagrange equations. 
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